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Axially  symmetric  experiments  produced  some  unexpected  results,  for 
example:  use  of  a  2  mm  jacket  of  lead  around  a  standard  pipe  produced 
holes  about  twice  as  great  in  depth  and  volume  as  the  average  of  the 
standard  configurations;  use  of  an  air-filled  pipe  with  wall  thickness 
of  0.7  mm  produced  triple  the  standard  hole  size.  Two-dimensional 
calculations  for  the  Standard  and  Lead  Wrap  cases  which  did  not  include 
the  thin  pipe  wall  were  totally  inconsistent  with  these  results. 

Calculations  including  the  wall  were  quite  consistent  with  the  data  if 
relative  energy  in  the  particulate  jet  of  wall  material  formed  on  collapse 
is  compared  with  measured  hole  volume.  In  terms  of  absolute  comparisons, 
the  calculated  results  are  consistent  only  if  the  cratering  efficiency 
of  the  shower  of  particulate  matter  is  assumed  to  be  three  to  five  times 
smaller  than  that  of  a  single  particle  at  the  same  velocity  -  which  is 
considered  reasonable. 

Simplified  calculation  including  the  wall  were  performed  for  all 
three  axisymmetric  configurations  in  the  frame  of  reference  of  a  steady 
shock  to  provide  guidance  for  the  calculations  of  the  actual  configurations 
with  a  decaying  shock.  These  produced  the  unexpected  result  that  all  three 
jet  particulate  mass  and  energy  at  equal  rates  as  a  steady  state  is 
approached.  This  result  is  explained  by  an  adaptation  of  impulsive  jetting 
theory.  The  steady  state  results  put  constraints  on  explanations  of  the 
laboratory  results.  The  stronger  jetting  in  the  Lead  Wrap  case,  compared 
with  the  Standard,  is  attributed  to  slower  response  of  the  region  of  the 
source  of  the  jet  to  decay  in  the  driving  pressure  in  the  wet  sand  -  the  Lead 
Wrap  configuration  is  better  tamped.  The  0.7  mm  air  filled  pipe  result  is 
explained  analogously. 

*  Configurations  employing  helical  spiral  asymmetries,  0.3  to  0.5  mm 
thick,  also  produced  unexpected  results.  Spirals  resting  against  the  outer 
pipe  wall  resulted  in  essentially  no  change  in  target  damage,  whereas 
spirals  resting  against  the  inner'wall  virtually  eliminated  target  damage. 

As  a  consequence,  original  intuition  that  the  spiral  would  prevent  cratering 
because  of  inertial  effects  was  abandoned.  A  postulated  mechanism  attribu- 
.  ting  the  spiral  effect  to  turning  of  the  unattached  spiral  into  the  pipe 
by  the  action  of  the  faster-arriving  vapor  flow  (generated  by  energetic  pipe 
implosion  near  the  source)  is  presented  in  detail.  This  mechanism  is 
shown  to  be  consistent  with  all  data  acquired  so  far  in  the  series  of  exper¬ 
iments. 
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SECTION  1 
INTRODUCTION 


A  series  of  high  explosive-driven  laboratory  experiments  has  been 
carried  out  by  Physics  International  Co.,  (PI)  over  the  past  several 
years  to  simulate  energy  injection  into  a  nuclear  1 ine-of-sight  pipe 
(LOS).  The  purpose  was  to  gain  better  understanding  of  the  mechanism 
of  energy  injection  into  the  LOS  and  to  test  means  of  minimizing  this 
energy  through  the  use  of  various  symmetric  and  asymmetric  configura¬ 
tions. 


The  experiments  considered  in  this  report  are  known  as  LS  II, 
III,  and  IV.  These  each  employed  a  spherical  charge  of  nitromethane 
embedded  in  a  large  container  of  saturated  sand.  This  drove  a  spher¬ 
ically  divergent  shock  in  the  sand  as  in  an  underground  nuclear  test. 

A  number  of  pipe  configurations  extending  in  various  directions  from 
the  sphere  were  employed  on  each  test  (nearly  all  were  evacuated). 

The  shock-driven  collapse  of  these  pipes  caused  energetic  jetting. 
Aluminum  witness  plates  at  the  end  of  each  pipe  were  used  to  assess 
the  jetting.  Figure  1.1  shows  the  arrangement  of  LS  II;  LS  III,  IV 
were  similar.  Greater  detail  is  available  in  Ref.  1. 

The  detailed  findings  of  these  experiments  are  described  in 
Refs.  1  and  2.  Ref  3  also  includes  an  adequate  presentation  of  LS  II, 
III  results.  Table  1.1  contains  a  brief  description  of  all  pipe 
configurations  tested  in  the  series  and  the  penetration  depth  and  hole 
volume  in  the  witness  plate  for  each.  The  major  findings  of  the  LS 
series  may  be  summarized  as  follows: 
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•  The  Standard  symmetric  configuration  was  employed  eight 
times  in  all.  A  relatively  wide  variation  in  target  damage 
was  observed. 

•  Use  of  an  external  spiral  asymmetry  or  an  internal  axial 
strip  asymmetry  produced  no  qualitative  reduction  in 
damage  to  the  target. 

•  Use  of  an  internal  spiral  asymmetry  drastically  reduced 
the  damage  to  the  target  for  all  spiral  materials  used 
(polyolefin,  steel,  lead). 

•  The  damage  reduction  remained  dramatic  for  spiral 
asymmetries  (or  for  spiral  or  ring  baffles  inside  the 
pipe)  occupying  only  the  last  38  cm  of  the  approximately 
1  m  total  pipe  length  but  was  only  moderate  for  spiral 
asymmetries  occupying  the  first  30  cm  or  less. 

•  Mufflers  occupying  the  last  38  cm  of  pipe  length  reduced 
the  average  damage  only  moderately,  whether  the  baffles 
employed  were  spirals  or  symmetric  rings  (these  baffles 
did  not  intrude  into  the  pipe). 

•  Use  of  a  heavy  symmetric  Lead  Wrap  with  the  Standard 
configurations  increased  the  average  damage  significantly 
-  considerable  variation  in  damage  was  again  observed  in 
duplicate  experiments. 

•  The  greatest  damage  level  of  all  was  attained  with  an 
air-filled  pipe  with  walls  2.3  times  as  thick  as  the 
Standard  pipe.  This  is  attributed  to  the  pipe  thickness  - 
a  repeat  of  the  Standard  with  air  produced  results  similar 
to  the  average  of  the  Standard  cases. 
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•  Air  at  1  atm  apparently  had  no  effect  on  the  performance 
of  any  of  the  symmetric  or  asymmetric  configurations  on 
which  it  was  used. 

•  "Standoff"  experiments  using  regions  of  wet  sand  of 
varying  length  between  the  spherical  source  and  the 
beginning  of  the  pipe  indicate  that  most  of  the  jetting 
energy  is  injected  in  the  first  25  cm  of  pipe  length  and 
nearly  all  of  it  in  the  first  50  cm. 

The  calculational  and  analytical  effort  by  S  Cubed  on  this 
program  has  two  principal  objectives:  first,  to  determine  how  well 
calculations  match  the  axially  symmetric  results  and  make  any 
necessary  improvements;  and  second,  to  provide  scenarios  to  explain 
both  the  symmetric  and  the  asymmetric  observed  phenomena.  The  initial 
effort  on  this  program,  which  was  only  partially  successful,  is 
described  in  Ref.  3.  The  effort  described  here  tentatively  fulfills 
the  program  objectives,  subject  to  the  findings  of  the  forthcoming  LS 
V  experiment  which  is  designed  In  part  to  clarify  the  proposed 
scenarios.  Nevertheless,  some  further  calculations  would  be  desirable 
as  time  and  resources  permit. 

The  calculations  and  analyses  of  the  symmetric  configurations 
are  discussed  in  Section  2.  A  scenario  for  the  performance  of  the 
asymmetric  configurations,  clarifying,  with  the  guidance  of  LS  IV 
data,  issues  that  were  left  unresolved  in  Ref.  3,  is  presented  in 
Section  3.  The  implications  of  the  results  for  underground  nuclear 
testing  are  discussed  in  Section  4.  A  summary  of  conclusions  and 
recommendations  is  given  in  Section  5. 
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SECTION  2 

SYMMETRIC  CONFIGURATIONS 


The  symmetric  configurations  used  in  the  experiments  produced 
some  unexpected  results  which  have  not  been  previously  explained. 

This  section  is  concerned  with  those,  in  particular,  the  comparisons 
between  the  Standard  configuration,  the  configuration  with  a  heavy 
Lead  Wrap,  and  the  "Pinex"  configuration  -  which  contained  air  in  the 
pipe  at  1  atm  (the  others  were  evacuated).  Results  which  were  already 
understood  in  terms  of  the  calculational  analysis  in  Ref.  3,  e.g.,  the 
"Standoff"  configurations,  are  not  considered  here. 

Calculations  of  these  conf igurations  indicate  that  shock-driven 
collapse  of  the  steel  pipe  causes  some  of  the  pipe  to  vaporize, 
resulting  in  an  energetic  flow  of  vapor.  These  also  indicate  the 
formation  of  a  jet  consisting  of  some  of  the  pipe  material  that  did 
not  vaporize.  The  volume  of  the  hole  in  the  witness  plate  is  con¬ 
sidered  a  good  comparative  indicator  of  the  energy  in  the  condensate 
jet  generated  by  the  pipe  collapse.  The  incident  vapor  flow  is  not 
thought  capable  of  contributing  much  to  the  excavation  process.  These 
statements  are  supported  later  in  this  Section.  Therefore,  calculated 
jet  energies  will  be  compared  with  measured  hole  volumes  as  a  test  of 
the  realism  of  the  calculations.  Table  2.1  lists  the  hole-volume 
statistics.  It  is  seen  that  a  considerable  band  of  outcomes  is 
observed  from  groups  of  seemingly  identical  experiments.  This  causes 
some  uncertainty  as  to  the  "real"  value  of  the  hole  volume,  especially 
for  the  Pinex  case,  of  which  only  two  were  employed  and  these  used 
quite  different  pipe  thicknesses. 

The  finding  that  the  jetting  in  the  Lead  Wrap  case  was 
considerably  more  energetic  than  the  Standard  was  quite  unexpected. 
This  configuration  was  thought  to  experience  slower  collapse  and  less- 
energetic  jetting.  However,  the  calculations  described  below  are 
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table  2.1  hole  volume  statistics  for  symmetric  configurations 


Standard 


Number  of  shots  8 

Pipe  Thickness,  cm  .030 


Average  hole 

volume,  cc.  17.8 

Standard 

deviation,  cc  7.8 

Experiment  numbers 

LS  II  -  1,2 
IS  III  -  1, 
10,11,14 
LS  IV  -  1,2 


_ P i nex _ 

Wr*P  Thick  Pipe  Thin  Pipe 


3 


1 


.030  .071 

+  .203  (lead) 


.030 


39.6  61.5  25# 


7.7 


LS  II  -  8  IS  Hi  _  20  LS  IV  -  17 
LS  IV  -  5,7 
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consistent  with  the  observations.  The  Pinex  (thick  pipe)  result  was 
even  more  surprising;  the  hole  volumes  for  Standard,  Lead  Wrap,  and 
Pinex  were  approximately  in  the  ratio  of  1  :  2  :  3.  However,  calcu¬ 
lations  and  approximate  theoretical  considerations,  described  in  Ref. 
3,  indicated  that  the  1  atm  air  present  in  the  Pinex  pipe  had  vir¬ 
tually  no  effect  on  the  outcome;  the  result  was  attributed  to  the 
greater  pipe  mass.  The  thin-pipe  Pinex  result  (LS  IV-17),  which  lies 
well  within  the  band  of  Standard  results,  bears  this  out.  Further 
discussion  of  the  mechanisms  responsible  for  the  various  results,  and 
their  implications,  will  follow  the  calculational  results. 

The  20  Eulerian  STREAK  calculations  described  in  Ref.  3  did  not 
include  the  thin  steel  pipe.  This  approximation  was  thought  to  be 
reasonable  at  the  time  and  allowed  substantial  savings  in  calcula¬ 
tional  costs.  The  findings  of  Ref.  3,  however,  suggested  that  the 
liner  (pipe)  must  be  represented.  Test  calculations  with  a  simple 
geometry  were  then  undertaken  to  determine  the  minimum  zonal  resolu¬ 
tion  required  to  give  an  adequate  description  of  liner  jetting. 

(These  calculations,  which  are  the  steady-state  results  described  in 
Appendix  A,  proved  so  informative  that  they  were  eventually  performed 
for  geometries  representative  Of  the  Standard,  Lead  Wrap,  and  Pinex 
configurations.)  It  was  found  that  the  zoning  requirements  were  less 
stringent  than  was  previously  thought.  This  is  evidently  a  conse¬ 
quence  of  the  use  of  a  dual  velocity  in  mixed  cells  in  recent  versions 
of  STREAK.  This  allows  slip  flow  along  material  boundaries  and 
effectively  minimizes  the  numerical  diffusion  responsible  for  an 
erroneous  slowing  of  vapor  adjacent  to  a  solid  wall  which  is  charac¬ 
teristic  of  Eulerian  codes.  It  also  has  the  beneficial  side  effect  of 
allowing  a  jet  to  initiate  without  being  5  or  6  zones  in  radial  extent 
-  this  limitation  had  been  experienced  with  Eulerian  codes  employing  a 
single  velocity  in  mixed  cells. 

The  Standard  and  Lead  Wrap  configurations  were  then  calculated 
with  the  initial  configurations  and  zoning  shown  in  Fig.  2.1.  Subse¬ 
quent  rezones  changed  the  axial  but  not  the  radial  zoning;  zoning 
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was  never  changed  in  the  immediate  vicinity  of  the  collapse  point  or  the 
ground  shock.  The  calculations  employed  the  "packaging"  feature  of 
STREAK  whereby  different  regions  of  the  problem,  requiring  different 
zonal  resolution,  are  calculated  in  separate  packages  which  communicate 
every  cycle.  In  this  manner  the  near-spherical  flow  some  distance  from 
the  pipe  could  be  treated  with  coarser  zoning  than  that  required  in  the 
pipe  itself. 

Figures  2.2  and  2.3  show  the  positions  of  the  tracer  strings 
defining  the  material  boundaries  at  various  times  in  the  calculations. 
3oth  calculations  indicate  the  formation  of  a  jet  of  condensed  material 
into  the  pipe  ahead  of  the  shock-driven  collapse.  Both  also  contain 
substantial  iron  vapor  in  the  pipe  which  was  generated  by  the  energetic 
collapse  of  the  pipe  wall.*  The  vapor  pressure  also  caused  pipe  expan¬ 
sion  ranging  to  about  15  in  the  Standard  calculation  (Fig.  2.2).  Very 
little  expansion  was  seen  in  the  Lead  Wrap  calculation  because  of  the 
great  mass  of  the  lead;  vapor  pressures  calculated  were  of  the  same 
order  as  those  in  the  Standard  case. 

It  is  seen  in  Figs.  2.2,  2.3  that  in  the  late  stages  of  each 
calculation  the  tracer  particles  defining  the  material  boundary  between 
the  steel  pipe  and  the  vapor  in  the  pipe  have  developed  a  distinct 
choppiness.  This  is  especially  true  at  128  usee  in  the  Standard  case, 
where  the  calculation  was  stopped.  The  source  of  this  difficulty  is 
that  the  thickness  of  the  undisturbed  pipe  is  only  about  one-third  of 
the  local  radial  zone  size.  Where  a  material  region  is  so  poorly 
resolved  the  algorithms  used  in  describing  tracer  velocity  and  mixed 
cell  pressure  will  occasionally  combine  to  produce  erratic  tracer  motion 


*  Liner  material  vaporized  by  collapse-induced  heating  is  transferred 
across  the  tracer  string  into  the  open  pipe.  Details  of  the  technique 
are  described  in  Ref.  9. 
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at  a  vapor/solid  boundary.*  After  examination  of  the  calculational 
edits,  we  believe  the  Standard  results  are  credible  until  just  before 
128  us  and  the  Lead  Wrap  results  until  at  least  199  us.  The  compari¬ 
son  of  jetting  energies  described  below  is  based  on  values  which 
leveled  off  before  these  times.  Furthermore,  the  condensate  jet  and 
the  vapor  are  formed  near  the  collapse  point,  away  from  the  tracer 
difficulty.  The  tracer  difficulty  involved  radial  velocities  of  at 
most  several  hundred  m/sec,  more  than  an  order  of  magnitude  lower  than 
the  axial  velocities  in  the  problem. 

Two  useful  quantities  for  comparison  are  the  calculated  energy 
in  the  condensate  jet  and  in  the  vapor.  These  are  shown  as  functions 
of  tir**-  in  Figs.  2.4  and  2.5.  The  plot  for  condensates  shows  kinetic 
energy  only;  their  internal  energy  is  much  smaller  and  is  not  thought 
to  contribute  to  penetration  of  the  witness  plate.  The  calculated  jet 
energy  *or  the  Lead  Wrap  levels  off  at  about  68.5  KJ  in  5.1  g,  vs  about 
..'•f  KJ  in  3.2  g  for  the  Standard.  The  vapor  energy  (internal  plus  kin¬ 
etic)  is  larger  for  the  Standard,  apparently  because  of  the  higher 
collapse  velocity.  The  calculated  jet  velocities  may  be  compared  with 
jet  velocities  inferred  from  the  data.  The  witness  plate  on  LS  1 1 1-14, 
a  Standard  configuration,  contained  diagnostic  pins  at  several  depths 
to  indicate  the  time  scale  of  growth  of  the  hole.  These  data  are 
shown  in  Fig.  2.6  along  with  vapor  arrival  data  from  this  and  several 
other  pipes  on  LS  III.  It  is  strongly  suggested  from  the  data  that 


*  These  algorithms  are  described  in  Ref.  12.  Briefly,  tracer  velo¬ 
cities  are  interpolated  from  cell  centers;  mixed  cell  pressure  is 
determined  by  volume-weighting  of  the  component  pressures.  Unless 
very  closely  spaced  tracers  are  used  to  describe  such  a  thin  strip  of 
material  as  the  pipe  wall,  radial  motion  may  at  times  cause  the  tracer 
strings  describing  opposite  sides  of  the  strip  to  attempt  to  cross. 

The  automatic  tracer  sorting  routine  responds  by  combining  the  opposite 
sides  locally.  This  results  in  a  high  density  in  that  material  which 
can,  in  extreme  cases,  produce  an  anomolously  high  cell  pressure  which 
"kicks"  some  of  the  solid  material  into  the  available  low-density 
region. 
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Figure  2.6  Arrival  and  penetration  data  from  Standard  and 
Polyolefin  spiral  configurations  of  LS  III. 
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the  hole  formation  is  too  late  to  be  attributed  to  the  vapor  arrival. 
It  is  made  by  something  which  travels  at  an  average  speed  of  ^bout  5 
Km/sec,  which  compares  well  with  the  calculated  velocity  of  about  5.5 
Km/sec  for  the  leading  edge  of  the  jet.  The  evidence  of  vapor  arrival 
in  the  Polyolefin  Spiral  configurations  LS  111-7,8  (at  about  the  same 
time  as  in  the  Standard  case.  Fig.  2.6),  which  produced  no_  significant 
damage  to  the  target  plate,  also  supports  the  claim  that  particulates 
make  the  hole.  The  calculated  vapor  front  arrival  is  also  indicated 
in  the  figure.  Although  the  calculation  was  not  continued  until  the 
formation  of  a  reflected  vapor  shock  at  the  witness  plate,  it  is 
evident  that  the  calculated  velocity  of  the  front  is  consistent  with 
the  arrival  data. 


It  is  postulated  that  the  condensate  jet  makes  the  observed  hole 
by  impact  cratering.  There  are  considerable  data  on  craters  made  by 
hypervelocity  impact  of  single  particles,  of  the  order  of  1mm  diameter 
or  larger  against  a  target.  In  the  present  experiments  the  impacting 
condensates  are  almost  surely  a  shower  of  much  smaller  particles, 
probably  both  solids  and  liquid  droplets.  Nevertheless,  it  is  reason¬ 
able  to  look  for  a  correlation  between  the  hole  volume  and  the  impact¬ 
ing  kinetic  energy.  The  calculated  and  experimentally-inferred  jetting 
velocities  lie  in  the  range  of  about  5  to  6  Km/sec.  Various  literature 
on  hypervelocity  impact  cratering^4-^  suggests  that  the  ratio  of 
hole  volume  to  impacting  energy  is  a  weak  function  of  particle  velo¬ 
city.  Therefore,  it  is  reasonable  to  require  the  ratio  of  jetting 
energies  calculated  for  two  configurations  to  be  comparable  to  the 
ratio  of  the  observed  hole  volumes.  The  average  hole  volumes  observed 


for  the  Lead  Wrap  case  are  about  2.2  times  the  average  observed  for 
the  Standard.  The  ratio  of  the  calculated  maxima  of  is  1.9, 
which  is  considered  to  be  well  within  the  experimental  uncertainty. 


A  quantitative  comparison  with  hypervelocity  data  is  also 
interesting.  For  iron  incident  on  aluminum,  it  may  be  inferred  from 
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empirical  relations  found  in  Ref.  4  that  the  ratio  of  hole  volume  to 
impacting  energy  is  V/E  =  2  to  3  cc/KJ.  This  would  lead  one  to  ixpect 
hole  volumes  of  about  70  to  110  cc  for  the  Standard  configuration  <snd 
140  to  200  cc  for  the  Lead  Wrap.  These  values  are  about  3.5  to  5  times 
the  observed  volumes.  This  is  not  surprising,  since  exact  correspon¬ 
dence  with  the  single  particle  data  is  not  to  be  expected;  linear 
superposition  is  not  credible.  If  the  calculations  are  taken  at  face 
value,  this  would  imply  that  a  shower  of  tiny  particles  has  cratering 
efficiency  several  times  lower  than  a  single  particle.  This  is  plausi¬ 
ble  in  view  of  interactions  between  the  particles  and  the  fact  that 
some  particles  will  collide  with  crater  ejecta.  The  principal  finding 
of  this  exercise  is  that  the  variation  in  calculated  jet  energy  for 
the  different  configurations  correlates  very  well  with  the  observed 
variation  in  target  damage.  The  variation  of  calculated  vapor  energy, 
on  the  other  hand,  does  not.  The  calculated  E  is  actually  lower 

Vdp 

for  the  Lead  Wrap  case.*  This  is  plausible  because  of  the  lower  col¬ 
lapse  velocity.  Vapor  energy  should  not  be  expected  to  correlate  well 
with  target  damage.  It  is  difficult  to  postulate  a  plausible  mechanism 
by  which  vapor  impact  would  make  a  hole  just  the  same  size  as  the  pipe, 
when  pipe  rupture  surely  accompanies  the  reflection  of  the  incident 
vapor  flow  at  the  witness  plate. 

The  calculational  results  are  seen  to  be  consistent  with  the 
experimental  results.  However,  it  is  necessary  to  try  to  understand 
the  results  in  terms  of  physical  mechanisms.  The  steady-shock  calcu¬ 
lations  for  the  configurations  (Appendix  A)  proved  quite  useful  in 
this  regard.  They  produced  the  unexpected  result  that  all  three 


*  The  original  expectation  that  the  Lead  Wrap  would  be  beneficial 
was  based  on  calculational  studies  indicating  that  the  use  of  extra 
mass  around  the  pipe  near  a  nuclear  source  would  result  in  lower-energy 
vapor  flow.  Symmetric  distribution  was  indicated  to  be  at  least  as 
good  as  asymmetric  distribution  (Refs.  13,14).  This  was  one  reason 
for  including  the  Lead  Wrap  configuration  in  the  experiments. 
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configurations  jet  energy  at  the  same  rate  when  driven  by  identica1 
steady  shock  conditions,  once  the  transient  features  have  decayed. 

This  result  was  then  explained  by  a  modification  of  Birkhoff's  simple 
jetting  theory  -  see  Appendix  A.  Therefore,  the  observed  differences 
when  driven  by  a  decaying  shock  must  be  due  to  nonsteady  effects. 
Furthermore,  any  explanation  must  be  consistent  with  the  steady  state 
findings.  As  a  result  of  this  requirement  an  explanation  that  origi¬ 
nally  seemed  pleasing  had  to  be  discarded  (the  scenario  of  Ref.  3  which 
resulted  from  an  inappropr iate  application  of  steady  state  jetting 
theory) . 

The  postulated  explanation  is  based  on  ground  shock  decay.  Decay 
is  rapid;  when  an  element  reaches  the  axis  the  driving  pressure  field 
is  significantly  lower  than  when  the  element  was  first  driven.  The 
pressure  experienced  by  the  element  at  convergence  on  axis  is  related 
as  much  to  the  initial  driving  conditions  as  to  the  current  pressure 
in  the  ground  around  the  element,  and  a  finite  time  is  required  to 
adjust  to  the  current  ground  pressures.  This  time  is  longer  in  the 
Lead  Wrap  than  in  the  Standard  because  of  the  thickness  of  the  imploded 
lead  and  the  lower  sound  speed  in  lead.  The  longer  duration  of  the 
high  convergence  pressure  in  the  Lead  Wrap  case  is  confirmed  in  the 
calculational  edits.  The  jet  formation  is  better  "tamped"  in  the  Lead 
Wrap  case  than  in  the  Standard. 

This  argument  also  explains  why  t* e  particulate  jet  is  more 
energetic  in  the  thick-walled  Pinex  (LSI  1 1-20)  than  in  the  Standard. 
However,  it  does  not  explain  why  the  jet  is  more  energetic  in  the 
Pinex  than  in  the  Lead  Wrap;  the  Pinex  would  appear  able  to  adjust  to 
reduced  driving  pressure  more  rapidly  than  the  Lead  Wrap.  Neverthe¬ 
less,  even  though  only  one  such  Pinex  configuration  was  used,  the  data 
suggest  a  high  probability  that  it  produces  more  energetic  jetting 
than  the  Lead  Wrap.  In  the  absence  of  a  Pinex  calculation  a  possible 
explanation  is  offered.  Because  of  the  smaller  implosion  velocity 
experienced  by  the  Lead  Wrap,  an  element  of  liner  initially  at  a  given 


station  reaches  the  axis  later  and  farther  behind  the  decayinn  shock 
front  than  the  corresponding  element  in  the  Pinex.  Consequently,  the 
average  driving  pressure  experienced  by  it  during  implosion  is  less 
and  the  jetting  energy  lower,  despite  the  longer  time  required  to 
reduce  the  convergent  pressure  to  that  of  the  surrounding  ground.  It 
would  be  very  useful  to  perform  a  calculation  of  the  Pinex  configura¬ 
tion,  when  time  and  budget  constraints  permit,  for  the  insight  it  is 
1 ikely  to  provide. 

In  conclusion,  the  calculations  and  phenomenological  explana¬ 
tions  now  appear  to  be  consistent  with  the  experiments.  Unfortunately, 
it  was  not  possible  to  obtain  such  satisfactory  agreement  without 
inclusion  of  the  thin  pipe  wall  in  the  calculations.  Implications  of 
this  for  the  nuclear  case  are  discussed  in  Section  4. 

Other  symmetric  configurations  are  worthy  of  comment  although 
they  have  not  been  calculated,  in  particular,  LSII-6,  9,  10  which 
employed  polyolefin,  glass,  and  cardboard  liners  inside  the  Standard 
steel  pipes  (Table  1.1).  These  liners  had  considerably  less  mass  than 
the  pipe  walls.  Therefore,  the  collapse  dynamics  should  be  very  simi¬ 
lar  to  that  of  the  Standard  case;  the  hole  volumes  in  the  witness 
plates  should  also  be  very  similar.  All  three  hole  volumes  lie  within 
the  spread  , f  the  Standard  results  (although  glass  and  polyolefin 
liners  produced  slightly  deeper  holes).  In  fact,  the  Standard  hole 
volume  statistics  available  at  the  writing  of  Ref.  3  (LSII,  III)  were 
virtually  indistinguishable  when  these  three  liner  cases  were  included: 
19.8  *  7.9  cc  excluding  them,  19.8  ±  8.2  cc  including  them.  These 
results  are  consistent  with  the  above  discussion.  It  is  not  clear,  of 
course,  what  role  the  liners  play  in  the  jetting  dynamics. 
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SECTION  3 

ASYMMETRIC  CONFIGURATIONS 

The  experiments  with  asymmetric  configurations  also  produced 
major  surprises.  The  failure  of  external  spiral  asymmetr ies  to  have 
any  significant  effect,  even  if  the  spiral  was  made  of  lead,  was 
completely  unexpected.  On  the  other  hand,  the  ability  of  a  thin  metal 
or  even  plastic  spiral  loosely  inserted  inside  the  pipe  to  suppress 
all  significant  damage  to  the  target  plate  was  equally  surprising. 

It  is  clear  that  the  internal  spiral  performs  its  "miracle"  by 
some  mechanism  quite  different  from  the  mass  asymmetry  which  was 
naively  thought  significant  in  the  early  days  of  this  research  effort. 
It  is  indicated  that  the  internal  spiral  prevents  jetting  by  either 
(a)  the  blockage  of  the  particulate  jet  at  its  source  in  the  pipe  col¬ 
lapse  region  or  (b)  the  scatter  of  jetted  particulates  by  the  spiral 
downstream.  The  latter  is  favored  by  the  body  of  data  summarized  in 
Table  1.1.  In  particular,  note  that  an  internal  spiral  in  the  last 
0.38  m  of  pipe  (LS  IV-8,9)  was  essentially  as  effective  as  those  going 
the  full  length,  whereas  spirals  limited  to  the  first  0.15  to  0.30  m 
were  only  partially  effective  .  "Standoff"  experiments  (IS  1 1 1-2 ,  3, 

4  and  LS  IV-20E)  suqgest  that  the  large  majority  of  the  energy  of 
particulate  jetting  is  injected  into  the  pipe  in  the  first  0.50  m  or 
less.  The  symmetric  calculations  suggest  the  same  thing.  Taken 
together,  these  facts  support  mechanism  (b). 

The  next  question  is,  how  can  such  thin  spirals,  resting  flat 
against  the  pipe  wall,  so  completely  protect  the  target  plate?  The 
difference  between  the  vapor  and  particulate  velocities  suggests  an 
answer.  Both  data  and  calculations  for  the  Standard  configuration 
indicate  velocities  of  about  10  Km/sec  for  the  vapor  front  and 
6  Km/sec  for  the  leading  edge  of  the  particulate  jet.  Thus,  the  vapor 
can  interact  with  the  spiral  for  a  time  before  the  particulates  arrive 
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and  this  time  increases  with  range  from  the  HE.  This  observation, 
combined  with  the  fact  (long  overlooked  in  the  attempt  to  understand 
the  observations)  that  the  spirals  are  loosely  inserted  in  the  pipe, 
riot  attached  to  it,  suggests  the  possibility  of  a  significant  effect 
of  vapor  on  the  spiral.  It  is  postulated  that  the  fast  vapor  stag¬ 
nates  against  the  leading  edge  of  the  spiral  and  forces  it  away  from 
the  pipe  wall.  This  then  becomes  an  obstruction  to  the  flow  which 
magnifies  the  twisting  tendency  of  the  aerodynamic  forces. 

It  is  easy  to  demonstrate  with  an  approximate  analysis  (see 
Appendix  B  for  details)  that  the  aerodynamic  forces  and  the  time 
available  are  sufficient  for  considerable  turning  to  occur  before  the 
arrival  of  the  slower  penetrating  materials.  The  vapor  arrival 
precedes  the  particulates  of  the  jet  at  the  end  of  the  pipe  by  about 
60uS  in  the  Standard  case.  Stagnation  pressures  up  to  about  10  kbars 
are  available  to  twist  the  spiral  insert.  It  is  estimated  that  between 
80  and  90  ps  is  the  time  required  to  turn  the  wide  face  of  the  spiral 
perpendicular  to  the  flow.  During  the  60  us  available  before  particu¬ 
late  arrival  this  analysis  suggests  the  spiral  will  turn  about  15°  and 
protrude  into  the  pipe  about  5  mm.  It  is  noteworthy  that  the  2  mm  x  2 
mm  lead  helical  baffles  in  the  last  0.38  m  of  pipe  in  test  IS  1 1 1— 12 
provided  almost  complete  protection  to  the  witness  plate. 

It  appears  that  the  relatively  thin  spiral  can  act  as  a  "meteor 
bumper"  to  deflect  or  vaporize  particulate  matter  which  would  penetrate 
the  target  plate  and  produce  a  noticeable  crater  if  earlier  impact  did 
not  occur.  It  is  believed  that  most  particulates  pass  near  the  wall 
at  some  stage  of  their  trajectories  to  the  witness  plate,  otherwise 
the  necessary  protrusions  would  need  to  be  a  larger  fraction  of  the 
pipe  radius. 

This  mechanism  is  consistent  with  the  observation  from  test 
IS  I V— 22  in  which  a  0.51  mm  steel  spiral  replaced  the  0.36  mm  plastic 
spiral  which  had  been  used  earlier.  This  experiment  used  a  spiral  in 
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^nly  the  first  0.30  m  of  pipe.  Data  extrapolations  from  plastic  spirals 
in  the  first  0.15  and  0.23  m  (LS  III-15,  16)  had  suggested  that  this 
limited  spiral  would  suppress  the  jet  completely.  However,  a  hole 
volume  of  nearly  half  the  average  produced  by  the  Standard  pipes  (and 
almost  as  large  as  the  smallest  hole  of  the  Standard  set)  was  actually 
observed  in  the  witness  plate.  It  is  suggested  that  this  result  fol¬ 
lows  from  the  relatively  high  density  of  iron  compared  to  plastic.  The 

analysis  of  Appendix  B  finds  that  the  rate  of  increase  of  turning 

-2  -1/2 

angle,  9,  is  proportional  to  At  oH  at  low  9,  where  a  is  the 
thickness  of  the  spiral  insert,  L  its  width  in  the  axial  direction, 
and  its  density.  As  a  result  the  steel  spiral  turns  through  a 
smaller  angle  in  the  limited  time  available.  The  smaller  turning 
probably  means  decreased  effectiveness. 

The  physical  blockage  of  the  pipe  by  the  warped  spiral  should 
have  some  effect  on  the  vapor,  in  addition  to  stopping  the  particu¬ 
lates.  Unpublished  data  from  LS  IV  (not  yet  fully  reduced  at  this 
time)  suggest  that,  while  first  arrival  of  the  vapor  at  the  target 
plate  is  only  slightly  delayed  by  the  spiral  (as  found  in  the  LS  III 
comparisons  cited  in  Section  2),  the  spiral  can  cause  a  sizeable 
reduction  in  the  vapor  impulse.  This  is  supported  by  data  from  the 
nuclear  event  Flora  (Section  4). 

In  summary,  the  proposed  scenario  explaining  the  role  of  the 
spiral  in  terms  of  an  aerodynamic  twisting  mechanism  appears  to  be 
consistent  with  the  entire  body  of  data  acquired  to  date  by  PI. 

The  planned  LS  V  experiment  will  contain  some  specific  tests  of  this 
scenario.  For  example,  there  will  be  pipes  in  which  the  spiral  it 
solidly  attached  to  the  wall.  Whether  the  ultimate  scenario  is  the 
one  proposed  or  some  other,  its  data  base  should  be  sufficient  to 
permit  more  confidence  in  it  than  is  now  possible. 
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SECTION  4 

IMPLICATIONS  FOR  UNDERGROUND  NUCLEAR  TESTING 


There  is  now  more  basis  for  conjecture  on  underground  test 
implications  than  existed  at  the  writing  of  Ref.  3  because  of  the 
execution  in  1980  of  the  nuclear  test  Flora  by  the  Los  Alamos  National 
Laboratory.  This  test  included  an  experiment  to  determine  whether  a 
spiral  insert  would  slow  the  fast  plasma  sufficiently  to  increase  the 
time  available  before  shock-induced  destruction  occurred  in  certain 
electro-optical  equipment.  Two  15  cm  diameter,  20  m  long  tubes  were 
installed.  These  were  identical  except  that  one  contained  a  steel 
spiral  insert  scaled  approximately  from  those  used  in  the  PI 
simulations. 

The  spiral  insert  had  a  significant,  favorable  effect  on  the 
plasma  flow  in  Flora  (Ref.  7).  Slifer  cables  showed  the  gradual 
development  of  a  lag  in  first  arrivals  with  range  in  the  asymmetric 
pipe.  Pressure  histories  recorded  at  18  m  indicated  plasma  arrival  at 
about  0.72  ms  in  the  smooth  pipe  with  a  peak  pressure  of  about  2.5 
kbars;  in  the  asymmetric  pipe  the  arrival  was  at  about  0.82  ms  with  a 
peak  of  about  1.3  kbars.  Only  the  first  plug*  was  penetrated  in  the 
asymmetric  pipe;  three  plugs  were  penetrated  in  the  smooth  pipe.  It 
had  been  suspected,  prior  to  the  shot,  that  the  spiral  would  expand 
with  the  pipe,  rupture  on  expansion,  and  not  be  as  effective  as  in  the 
PI  configurations.  This  was  apparently  not  the  case.  The  key  to  the 
success  of  the  spiral  on  Flora,  as  in  the  PI  cases,  may  well  be  that 
the  spiral  was  not  attached  to  the  pipe. 


*  At  the  end  of  each  pipe  was  a  sequence  of  three  plugs  (A1  plates) 
at  about  23  cm  intervals,  the  first  1.27  cm  thick,  the  second  2.54  cm, 
and  the  third  3.81  cm.  Penetration  was  indicated  by  shorting  of  the 
TDR  cable  behind  the  plug  which  could  be  caused  by  vapor  pressure  or 
by  failure  of  the  plug. 
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The  primary  observed  effect  of  the  Flora  spiral  was  on  plasma 
flow.  Its  effect  on  particulate  jetting  is  unknown,  and  the  contribu 
tion  of  particulates  to  plug  penetration  is  not  known.  On  typical  HLOS 
tests  particulate  jetting  is  thought  to  be  an  even  more  severe  threat 
to  the  samples  than  the  plasma.  Therefore,  it  is  important  that 
nothing  impede  the  spiral's  function  as  a  stopper  of  particulates. 
According  to  calculations,  the  energetic  plasma  flow  in  the  LOS  causes 
substantial  pipe  expansion  while  the  flow  in  the  PI  simulations  does 
not.  Calculated  LOS  expansion  in  the  nuclear  case  is  usually  well  in 
excess  of  the  30%  believed  to  cause  pipe  rupture.  A  concern  with  the 
use  of  spirals  in  LOS  tests  is  that  the  ability  of  the  spiral  to  stop 
particulate  jetting  has  not  been  demonstrated  when  major  pipe  expan¬ 
sion  occurs.  Therefore,  the  normal  requirement  to  minimize  energetic 
plasma  flow  appears  completely  compatible  with  the  use  of  a  spiral. 
Obviously,  there  is  a  point  at  which  further  reduction  in  the  plasma 
flow  would  be  undesirable,  if  the  spiral  works  as  we  suggest.  But 
there  is  no  evidence  that  this  level  is  attainable  in  practice.  The 
prompt  plasma  flow  for  Diablo  Hawk  was  the  most  benign  ever  calculated 
or  measured,  but  the  calculations  indicate  that  it  caused  appreciable 
pipe  expansion  while  the  PI  simulations  did  not.  Specific  applications 
of  the  approximate  "twisting"  model  of  Appendix  B  suggest  that  the 
prompt  plasma  flow  calculated  for  Diablo  Hawk  could  turn  a  spiral 
perpendicular  to  the  flow  prior  to  the  arrival  of  particulate  jetting, 
at  least  to  a  range  of  the  order  of  50  m.  (This  assumes  a  spiral 
whose  width  equals  its  radius.) 

The  symmetric  Lead  Wrap  results  were  quite  unexpected  because  of 
earlier  classified  parametric  studies  indicating  the  superior  perfor¬ 
mance  of  dense,  long  extensions.  The  earlier  work  used  plasma  energy 
as  the  primary  figure  of  merit;  the  threat  associated  with  the  Lead 
Wrap  configuration  is  thought  to  be  jetting  particulates.  (Recall 
that  the  Lead  Wrap  calculation  described  in  Section  2  predicts  vapor 
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energy  30a'  less  than  the  Standard  calculation  as  well  as  later  arri¬ 
vals,  which  is  consistent  with  preliminary  pressure  history  data  from 

(2) 

the  LS  IV— 5.7  Lead  Wrap  pipes  . )  The  question  arises  whether  the 
Lead  Wrap  results  have  adverse  implications  for  dense  extensions.  The 
answer  appears  to  be  that  they  do  not.  Even  the  longest  extensions 
used  to  date  are  limited  to  the  near-source  region,  whereas  the  PI 
lead  wrap  went  the  full  length  of  the  pipe.  A  well-designed  extension 
is  overtaken  by  the  ground  shock-driven  closure  while  the  energy  in 
the  LOS  is  still  low.  The  ground  shock  pressures  at  the  ends  of  the 
Mighty  Epic,  Diablo  Hawk  and  Miners  Iron  tubal loy  extensions  were 
approximately  the  same  as  the  driving  pressure  at  the  PI  source.  Hence, 
the  two  situations  are  not  directly  comparable.  The  calculation-based 
understanding  of  the  LOS  results  is  consistent  with  data,  as  is  the 
calculation-based  understanding  of  the  PI  symmetric  results.  Most 
importantly,  dense  extensions  have  performed  well.  Therefore,  the 
conclusions  of  the  earlier  parametric  study  are  still  thought  to  be 
valid.  It  is  speculated  here  that  a  short  Lead  Wrap  near  the  source 
in  a  PI  simulation  would  result  in  no  more  jetting  than  the  Standard 
case  and  quite  possibly  less. 

Certainly,  extension  designs  have  not  yet  been  optimized.  The 
Miners  Iron  LOS  flow  was  more  energetic  than  expected.  It  is  recom¬ 
mended  that  further  optimization  studies  include  some  LOS  closure 
calculations  which  represent  the  thin  pipe  wall  which  appears  to  be 
important  in  the  PI  simulations.  This  would  be  useful  in  determining 
whether  the  use  of  dense  extensions  to  minimize  plasma  flow  causes 
locally-enhanced  particulate  jetting.  If  it  does,  a  spiral  insert 
should  be  able  to  control  it.  It  is  not  clear  that  greater  realism 
would  be  achieved  by  including  the  pipe  wall  beyond  the  extension 
unless  pipe  expansion  is  insufficient  to  rupture  the  pipe.  In  a  case 
where  major  rupture  is  indicated,  the  particulate  jet  could  in  reality 
consist  mostly  of  the  exposed  grout. 
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SECTION  5 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  are  based  on  the  experimental  data  and 
insight  gained  from  the  calculations: 

•  The  damage  to  the  witness  plate  appears  to  be  due  to  the 
jetting  particulate  matter,  the  primary  source  of  which  is 
the  steel  pipe. 

•  2D  Eulerian  calculations  of  the  symmetric  configurations 
appear  to  produce  jetting  results  consistent  with  the 
comparative  damage  levels  observed  on  the  Standard  and 
Lead  Wrap  experiments,  provided  the  thin  steel  pipe  wall 
is  explicitly  included  in  the  calculation. 

•  The  success  of  internal  spiral  asymmetries  in  preventing 
damage  to  the  witness  plate  and  the  relative  damage  levels 
observed  in  various  configurations  using  such  asymmetries 
in  only  a  portion  of  the  pipe  are  consistent  with  the 
postulate  that  the  spiral  is  separated  from  the  wall  by 
the  precursor  vapor  flow  and  acts  as  a  "meteor  bumper"  to 
deflect  or  vaporize  hypervelocity  particulate  matter. 

•  The  use  of  internal  spiral  asymmetries  in  nuclear  LOS 
events  has  great  potential  for  drastically  reducing  damage 
to  closures  or  experimental  objects  caused  by  particulates 
and  for  eliminating  the  need  for  other  hardware  to  prevent 
such  damage.  It  also  has  significant  potential  for 
delaying  and  reducing  plasma  flow. 

In  addition  to  the  above,  the  following  observation  has  been  so 
persistently  valid  throughout  this  work  that  it  warrants  inclusion  in 
the  above  list: 
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•  The  understanding  of  experimental  results  is  greatly 
enhanced  by  the  attempt  to  perform  calculations,  even  if 
these  are  initially  inconsistent  with  the  data;  and  the 
results  of  calculations  of  complicated  phenomena  are 
highly  suspect  unless  they  have  been  compared  with 
relevant  data.  Experimental  and  theoretical  work  advance 
hand  in  hand.  The  greatest  progress  comes  in  attempting 
to  reconcile  experiments  and  calculations;  for  example, 
the  spiral  mechanism  now  appears  to  be  very  different  from 
that  initially  imagined  when  inertial  effects  were  expected 
to  dominate. 

The  following  recommendations  are  also  offered: 

•  An  internal  spiral  asymmetry  should  be  tried  on  a  nuclear 
LOS  event  as  soon  as  possible.  Its  design  should  be 
carefully  based  on  the  experimental  findings  and 
phenomenological  understanding  obtained  from  the  PI 
simulations. 

•  Continuing  calcul ational  studies  of  nuclear  LOS  closure 
performance  should  include  some  calculations  in  which  the 
thin  pipe  wall  is  explicitly  represented. 
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APPENDIX  A 


CALCULATIONS  OF  JETTING  WITH  A  STEADY  SHOCK 


Calculations  of  jet  formation  in  the  frame  of  reference  of  an 
assumed  steady  shock  wave  were  initially  undertaken  as  an  economical 
way  to  determine  the  zoning  requirements  for  calculations  of  the  actual 
decaying-shock  configuration  including  adequate  resolution  of  the 
liner.  The  results  proved  to  be  so  useful  that  the  Standard,  Lead 
Wrap,  and  Pinex  configurations  were  all  calculated  in  this  manner. 

These  calculations,  performed  with  the  STREAK  code,  produced  some 
unexpected  findings,  which  were  later  explained  in  terms  of  a  simple 
theoretical  model.  The  model  even  gives  reasonable  quantitative  agree¬ 
ment  with  the  calculations.  Any  scenario  developed  for  the  performance 
of  the  actual  decaying-shock  configurations  must  be  consistent  with  the 
insights  derived  from  this  exercise. 

Conditions  selected  for  the  calculation  were  taken  from  a 
previous  calculation  for  the  PI  Standard  configuration  just  beyond  the 
range  where  pipe  collapse  on  axis  first  occurred.  (That  calculation 
did  not  include  the  steel  liner.)  The  shock  conditions  were: 
pressure  P  =  119  kbar,  density  o  =  2.7187  g/cc  (wet  sand), 

S  S  g 

specific  internal  energy  e$  =  7.9  x  10  erg/g.  With  ambient  den¬ 
sity  o0  =  2.0198  g/cc  this  corresponds  to  a  shock  velocity  of 
V$  =  4.785  km/sec.  The  material  behind  the  shock  has  velocity 
u'=  3.555  km/sec  in  the  shock  frame  of  reference.  Fig.  A.l  shows  the 
configuration  and  the  zoning.  The  above  values  of  o$,  e$,  and  u' 
were  imposed  at  the  top  and  right  of  the  grid  at  all  times.  In  this 
frame  of  reference  and  orientation  the  jet  goes  out  the  top  of  the 
grid. 
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PrtECJsDli<G  PaOE  BUjNK-WOT  FILLED 


The  zoning  was  employed  on  a  trial  basis  and  was  not  expected  to 

be  adequate  for  jet  formation.  Previous  experience  of  other  invest.iga 

tors  with  jetting  calculations  using  the  HELP  code  had  produced  a  -ule 

of  thumb  requiring  at  least  5  or  6  zones  radially  across  the  jet  if  it 
(8) 

is  to  form.  ' .  The  present  finding  that  two  zones  is  sufficient 
has  tentatively  been  attributed  to  the  use  of  a  "slip  velocity"  in 
mixed-material  cells  in  STREAK;  each  material  has  its  separate  velo¬ 
city.  This  minimizes  the  numerical  diffusion  inherent  in  Eulerian 
codes. 

Pour  calculations  were  carried  out:  the  Standard  configuration, 
the  Pinex  conf iguration  (thicker  wall,  plus  1  atm  ambient  air  in  the 
pipe),  the  Lead  Wrap  configuration,  and  a  more-finely  zoned  version  of 
the  Lead  Wrap  -  to  assess  the  influence  of  zoning.  The  calculations 
were  started  with  the  pipe  at  its  original  radius  everywhere  and 
carried  until  a  steady  flow  configuration  with  jetting  was  achieved. 

It  is  useful  to  compare  the  mass  and  energy  moving  ahead  of  the 
collapse  point  in  these  calculations.  After  steady  flow  is  achieved, 
the  rates  of  mass  and  energy  flux  out  the  top  of  the  grid  may  be 
compared.  (This  is  a  more  convenient  quantity  to  edit).  Figs.  A. 2 
and  A. 3  show  the  mass  in  the  jet  and  the  mass  in  the  vapor  that  have 
gone  out  the  top  as  functions  of  time.  The  vapor  is  made  at  the 
collapse  point.  Liner  material  is  vaporized  as  required  by  collapse- 
induced  heating.  The  mass  vaporized  in  this  manner  at  each  cycle  is 
transferred  across  the  tracer  particle  string  defining  the  material 
boundary  (technique  described  in  Ref.  9).  It  is  noteworthy  that  the 
asymptotic  mass  flux  rate  into  the  jet,  M. ,  is  nearly  the  same  in 

J 

all  calculations.  This  result,  which  was  at  first  surprising,  is 
consistent  with  a  simplified  theory  of  jetting  -  see  below.  (The  fine 
zoned  Lead  Wrap  case  was  not  carried  far  enough  for  this  comparison; 
it  was  useful  only  in  demonstrating  that  the  coarser  zoning  provided 
an  adequate  description  of  the  jetting). 
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The  plots  of  energy  out  the  top  of  the  grid  in  the  jet  and  in 
the  vapor.  Figs.  A. 4  and  A. 5,  have  the  same  character  as  the  mass  flux 
plots.  The  energy  flux  rate,  E-,  for  the  condensate  approaches 

J 

nearly  the  same  value  in  all,  again  consistent  with  the  simplified 
theory. 

In  order  to  explain  the  phenomenon  of  equal  M.  and  E.  for 

J  J 

different  configurations,  it  is  necessary  to  digress  and  review  the 
impulsive-jetting  theory  of  Birkhoff,  et  al^’^.  Refer  to  Fig. 

A. 6  which  depicts  a  detonation  wave  propagating  through  a  high  explo¬ 
sive  surrounding  a  thin  liner.  A  steady  state  exists  in  the  frame  of 
reference  moving  with  the  apex.  In  this  frame  u£  is  the  detonation 
velocity,  the  velocity  of  material  toward  the  apex,  and  vq  the 
impulsive  velocity  imparted  to  the  liner  by  the  HE.  The  "impulsive" 
approximation  is  based  on  the  rapid  relief  of  pressure  behind  the 
detonation  wave  because  of  expansion,  v  is  the  velocity  of  the 

a 

apex  in  the  laboratory  frame.  For  simplicity,  we  need  consider  only 
the  cylindrical  pipe  case  (a  *  0). 

The  following  expressions  apply  (see  vector  diagram  in  Fig.  A. 6) 


vq  =  2u£  sin  ? 


V1  " 


(A.  1) 
(A.  2) 


va  -  uJl 


(A.3) 


The  steady  state  assumption  permits  the  Bernoulli  Equation  to  be  used: 

.2 


P  u 

e  +  -  +  ■*-  =  constant 

o  i 


(A. 4) 
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Figure  A. 4  Flux  of  condensate  kinetic  energy  out  of 
the  computational  grid. 
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Figure  A. 5  Flux  of  total  vapor  energy  out  of  the 
computational  grid. 
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Figure  A. 6  Impulsive  jetting  configuration  considered 
in  simplified  theory. 
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Since  e  is  changed  but  little  by  the  impulsive  push  and  P  =  0  on  any 
free  surface,  all  free  surface  velocities  are  approximately  u«  in  the 
apex  frame  of  reference.  This  applies  to  the  jet  and  to  the  rearward- 
moving  slug.  Therefore,  in  the  shock  frame,  the  slug  and  jet 
velocities  are 


vs  -  -u£  (A. 5) 

Vj  =  u£  (A. 6) 

In  the  laboratory  frame,  vs  =  0  and  vj  =  2u£- 

Straightforward  application  of  conservation  of  mass  and  axial 
momentum^  gives  the  following: 

2  a 

mj  =  m  sin  |  (A. 7) 

2  a 

m$  =  m  cos  j  (A. 8) 


where  m  is  the  liner  mass  per  unit  axial  length  and  m.,  m  are  the 

J  — 

components  of  m  entering  the  jet  and  the  slug.  The  energy  fed  into 
the  jet,  per  unit  axial  length  of  pipe,  is: 


e 


j 


(A. 9) 


(m  vq/2  in  the  laboratory  frame). 

The  above  results  for  impulsively-driven  liner  jetting  may  be 
adapted  to  the  steady  shock-driven  liner  problem.  This  is  different 
in  that  the  pressure  field  responsible  for  collapse  remains  substan¬ 
tial  well  behind  the  shock  and  the  liner  continues  to  accelerate  until 
it  reaches  the  axis  (Fig.  A. 7).  The  collapsing  liner  is  not  conical. 
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Figure  A.  7  Steady  shock-driven  jetting  configuration. 
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It  is  useful  to  define  s  as  the  total  angle  turned  by  the  streamline 
into  the  stagnation  point  (Fig.  A. 7).  The  stagnation  point,  whicn  is 
stationary  in  the  moving  frame  of  reference,  now  plays  the  role  of  the 
apex.  The  Bernoulli  Equation  still  holds  for  the  liner.  We  again 
obtain  the  result  that  all  free  surface  velocities  have  the  magnitude 
u^.  (Actually,  there  is  substantial  pressure  in  the  pipe,  but  after  a 
steady  jet  is  established  this  pressure  is  fairly  uniform  between  the 
jet  and  the  shock  position,  and  changes  of  e  have  been  found  to  be 
unimportant,  so  the  result  still  holds.)  Eq.  A.l  still  holds  if  v’ 
is  the  total  velocity  change  imparted  by  the  post-shock  pressure  field 
and  if  the  angle  of  the  stagnation  streamline  does  not  change  much  in 
the  region  around  the  stagnation  point  where  pressure  changes  are 
appreciable.  It  follows  that  the  entire  set  of  Eqs.  A.l  -  A. 9  holds. 

Equation  (A. 9)  is  helpful  in  understanding  the  result  of  nearly- 
equal  E.'s  for  the  various  configurations.  It  states  that  the  jet- 

J 

ting  energy  is  equal  to  the  kinetic  energy  of  collapse  imparted  by  the 
driving  forces.  In  the  steady  shock  case  this  should  be  proportional 
to  the  driving  pressure  times  the  collapse  distance,  P$r0»  which 
is  the  same  in  all  calculated  configurations. 

Now  the  results  of  nearly-equal  M.'s  can  be  understood.  From 

J 

Eqs.  A. 6  and  A. 9: 


or 


2ej  2Ej 

Mj  “  mjV  =  ^2 

J  J  *  u*  ug, 


(A. 10) 


Since  the  E.'s  and  u5  's  are  the  same  for  all,  so  must  be  the  M.'s. 

.1  ^  J 
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Now  the  simple  theory  can  be  compared  with  the  calcul atior.al 
results.  Typical  results  from  each  calculation  after  a  quasi-steady 
flow  has  evolved  are  shown  in  Figs.  A. 8  and  A. 19.  It  should  first  be 
pointed  out  that  in  the  calculations  the  wet  sand  stemming  and  the 
liner  were,  for  simplicity,  fed  into  the  grid  at  the  same  velocity 
=  u/ =  3.555  km/sec,  the  velocity  of  the  shocked  sand  in  the  shod 
frame  of  reference.  This  would  not  occur  in  a  real  configuration. 


First  we  verify  that  the  solutions  satisfy  Eq.  A. 6.  The 
velocity  vector  plots  for  the  three  basic  calculations  are  shown  in 
Figs.  A. 11,  A. 15  and  A. 19.  On  each  of  these  is  shown  an  arrow  corres¬ 
ponding  to  an  axial  velocity  u^  .  Eq.  A. 6,  in  the  moving  frame  of 
reference,  requires  that  v^.  =  u^ .  It  is  evident  that  the  flow  in 
the  region  of  the  jetting  liner  most  obvious  in  the  tracer  plots 
(location  given  in  Figs.  A. 8,  A. 12,  and  A. 16)  satisfies  this 
relationship  to  reasonable  accuracy. 


It  is  not  possible  with  the  simplified  theory  to  make  an  a  priori 
prediction  of  the  rates  of  mass  and  energy  flow  associated  with  jet¬ 
ting.  There  is  not  sufficient  information  to  determine  B.  However, 
various  checks  on  the  internal  consistency  of  the  theory  can  be  made. 
The  rates  of  mass  and  energy  jetting  compared  with  the  theory  will  be 
the  totals  for  the  vapor  plus  the  particulate  jet;  the  theory  con¬ 
siders  a  single  jet  velocity  with  no  vapor,  whereas  the  calculation 
allows  different  velocities  for  the  vapor  and  the  condensed  phase. 
First,  values  of  8  will  be  inferred  from  calculated  values  of  M.  and 

J 

compared  with  values  estimated  from  the  tracer  and  contour  plots. 

Table  A.l  shows  the  relevant  quantities.  Given  m  and  the  calculated 

M.,  Eq.  A. 10  implies  a  value  of  m.  and  Eq.  A. 7  implies  a  value  of 
J  J 

B.  This  B  is  seen  to  compare  reasonably  well  with  values  inferred 
from  the  tracer  plots.  The  difference  is  the  greatest  for  the  Lead 
Wrap  case,  where  the  thin  liner  approximation  is  most  suspect.  In 
addition,  the  asymptotic  values  of  M,  estimated  from  Figs.  A. 2  and 

J 

A. 3,  and  especially  the  values  of  s  estimated  from  the  contour  and 
tracer  plots,  are  only  approximate. 
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Table  A.l  Comparison  of  collapse  angles, 8°,  implied  by  simplified 
_ theory  and  measured  from  calculation  output _  _ 


1 

1 

m(g/cm) 

j  Calculated  M, (kg/sec) 

mj(g/cm)jlnferred  6° 

6° Estimated  from 

Condensates 

.... 

-  Mj/u^  Eq.  A. 7 

Contours./T  racers 

1 

jstandard 

1.413 

212 

91 

0.852  102 

■  i 

| 

90 

1 

'Lead  Wrap 

16.68 

212 

37 

249 

0.700  ]  24  j 

31 

Pinex 

i 

3.225  j 
1  1 

212 

37 

- 1 

249 

0.7C0  56 

i  ! 

58 

A  value  of  £.  can  be  inferred  from  a  calculated  M.  using  Eq. 

J  »  J 

A. 10.  This  can  be  compared  with  the  E.  edited  from  the  calcula- 

0 

tions.  The  comparisons  are  shown  in  Table  A. 2.  The  inferred  values 
support  the  qualitative  validity  of  the  simplified  theory  although 
they  are  25%  to  28%  low.  The  effect  of  partial  vaporization  of  jetted 
material  on  jetting  efficiency  is,  of  course,  not  addressed  by  the 
theory.  It  is  not  unreasonable  to  expect  that  this  would  allow  more 
energy  to  jet  because  vapor  can  escape  the  stagnation  region  into  the 
open  pipe  around  the  narrow  condensate  jet.  The  calculated  result  is 
consistent  with  this  expectation. 


Table  A. 2  Comparison  of  jetting  energy  rate  inferred  from 
_ Mj  and  calculated  directly _ 


1 - 

| 

Directly  Calculated 

Ei 

1 

Inferred  Ej(GJ/sec) 

1.91 

Condensates 

Vapor 

Total 

i 

1  Standard 

1.59 

1.08 

2.67 

Lead  Wrap 

i 

i 

1.57 

1.59 

0.50 

i 

2.09 

Pinex 


1.57 


1.59 


0.50 


2.09 


Axial  DisLance,  cm 


Figure  A. 8  Tracer  positions  in  idealized  Standard 
configuration,  25  ns. 
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Axial  Distance,  cm 


i 


Radial  Distance,  cm. 


Figure  A. 9  Pressure  contours  in  idealized  Standard 
configuration,  25  fj. s. 
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Figure  A. 11  Velocity  vectors  in  idealized  Standard 
configuration,  25  ps. 
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Figure  A. 15  Velocity  vectors  in  idealized  Lead  Wrap 
Configuration,  25  jis . 
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Figure  A. 16  Tracer  positions  in  idealized  Pinex 
configuration,  25  jus. 
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Figure  A. 17  pressure  contours  in  idealized  Pinex 
configuration,  25  /is. 
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Figure  A. 19  Velocity  vectors  in  idealized  Pinex 
configuration,  25  ps . 
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Eq.  A. 4  can  be  used  to  estimate  the  pressure  at  the  stagnation 
point.  Ignoring  the  small  pressure  and  internal  energy  in  the  liner 
as  it  feeds  into  the  jet,  and  denoting  stagnation  conditions  with 
subscript  s,  the  stagnation  pressure  is  approximately 


Ps  = 


ps( 


~r 


-  e. 


(A. 11) 


The  stagnation  pressure  occurs  in  the  collapsed  iron  liner.  The 
following  simple  equation  of  state  for  iron  may  be  used: 


P  =  Ay  +  By2  ,  v  =  ■£-  -1 
°o 


(A. 12) 


where  A  =  1.28  Mbar,  B  =  1.05  Mbar,  oQ  =  7.87  g/cc.  e^,  o$,  and 
P$  can  be  found  by  a  trial  and  error  method  satisfying  Eqs.  A. 11, 

A. 12,  and  e$  = Jpd{llo)  assuming  an  adiabatic  compression  to  the 
stagnation  state.  For  u^  =  3.555  km/sec  the  result  is  Pg  *  588  kbar. 
The  calculations  all  produce  results  fluctuating  between  about  350  and 
700  kbar,  which  is  typical  of  solutions  with  Eulerian  codes  at  a  point 
convergence  on  axis.  The  agreement  between  the  calculations  and 
the  estimate  with  Eqs.  A. 11,  A. 12  is  considered  good. 


It  is  concluded  that  the  simple  theory  is  quite  consistent  with 
the  calculational  results,  despite  its  simplicity.  In  particular,  it 
does  not  consider  the  effect  of  partial  vaporization  of  the  liner,  and 
it  assumes  a  thin  liner.  The  latter  assumption  is  highly  questionable 
in  the  Lead  Wrap  case.  This  is  also  the  case  where  the  inferred  and 
calculated  s's  are  in  greatest  disagreement  (about  25%).  It  is  impor¬ 
tant  to  emphasize  that  only  qualitative  comparisons  with  simple  theory 
are  possible  because  the  asymptotic  values  of  the  calculated  M.  and 
Ej  are  known  only  approximately  and  there  is  considerable  uncertainty 
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in  the  value  of  8°  estimated  with  the  aid  of  tracer  and  contour 
plots  -  especially  in  the  Lead  Wrap  case.  The  former  difficulty  can 
be  reduced  by  running  the  calculations  longer;  the  latter  cannot. 


It  is  also  concluded  that  the  zoning  employed  in  these  steady 
shock  frame  calculations  is  adequate  for  the  calculation  of  jetting 
with  a  spherically  diverging  shock. 


62 


APPENDIX  B 


AERODYNAMIC  TWISTING  OF  THE  HELICAL  SPIRAL 

The  stagnation  of  the  flow  against  the  forward  face  of  the 
spiral  generates  a  torque  about  the  point  i?i  Figure  B.l.  This 
causes  an  angular  displacement,  0.  As  9  increases,  the  torque 
increases  rapidly  and  so  does  the  angular  acceleration,  9. 

The  following  analysis  provides  an  estimate  of  the 
characteristic  time  required  for  9  to  become  very  large.  Rigid  body 
rotation  about  the  point  0  is  assumed.  This  tends  to  underestimate 
the  rate  of  turning  because  the  spiral  seems  likely  to  bend  about  some 
intermediate  point  along  the  length  L  and,  therefore,  to  have  in 
effect  a  smaller  moment  of  inertia.  On  the  other  hand,  the  omission 
of  strength  effects  tends  to  overestimate  the  turning  rate.  We  feel 
that,  on  balance,  the  estimate  is  pessimistic  (an  underestimate) 
because  the  strength  effect  is  small;  this  will  be  justified  later. 
Slippage  of  the  point  0  downstream  along  the  wall  is  also  likely  to 
occur  but  this  should  not  change  the  conclusions. 

The  moment  of  inertia  (I)  and  the  torque  (T)  per  unit  depth  into 
the  diagram  of  Figure  B.l  may  be  expressed  as  follows: 

.,3 
phaL 


2 

T  =  [L  sin  9  +  a  cos  9]2 


(B.l) 

(B.2) 


The  angular  acceleration  is 
9  =  T/I 


(B.  3) 
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(a)  Initial  configuration 


{b)  After  start  of  twisting 


Figure  B.l  Geometry  of  twisting  spiral 


Since  this  analysis  has  only  qualitative  validity,  it  is  justifiable  to 
simplify  the  problem  by  linearizing  in  0.  This  will  give  a  readily 
integrable  equation  for  0(t)  which  is  valid  for  small  0.  The  solution 
diverges  at  time  t„  ;  t*,  will  be  taken  as  an  estimate  of  the 

characteristic  time  for  0  to  become  ''large". 


g-  =  3  OU_ 


7  — ~7  [L  0  +  Al< 

°HaL 

For  convenience,  let 


a  _  3  pu‘ 

~  7  p^aL 


Then 


K  A 
b  =  j- 


9  =  a(0  +  b/ 

With  the  substitutions 

d  =  8  +  b,  t  =  ,/a"  t 

t".  42 


(8.4) 


The  initial  conditions  are 

4  (0)  =  0,  d(0)  =  b. 

The  initial  condition  on  4  will  be  satisfied  only  approximately. 


Two  solutions  are  available: 


(8  *t)‘ 


*  =  ? 


2A 


(8  *t  )• 


4'  - 


6A 


(8  ±T 


where  A  and  B  are  constants.  The  difficulty  is  that  these  cannot  be 
superposed  to  satisfy  the  initial  conditions.  The  differential 


65 


equation  requires  that  A  =  6.  Assuming  B>0,  the  lower  sign  must  be 
used  so  that  <b' >  0  as  required  physically.  <t>‘  (0)  >  0,  in 
contradiction  to  the  requirement,  but  it  is  sufficiently  small  to  be 
acceptable  (this  will  be  verified  later).  From  the  other  initial 
condition. 


*(0)  =  b  =  ^ 

It  follows  that 

B  =  (6L/a)1/2  , 

*'(0)  =  (§)1/2  (£)3/2 


9(0)  =  (*-)1/2  ^ 
°H  L2 


0  = - ^ ,  -  a/L 

(1  -  t/tj2 

where 

oo  u  '  0 


(B.5) 

(B.6) 

(B.7) 


0  >oo  as  t  >  t_  . 

oo 

Equations  (3.5)  and  (B.7)  yield  the  following 
0(O)t^  =  2a /L 

i.e.,  if  the  unwanted  initial  angular  velocity  9(0)  inherent  in  the 
above  solution  were  sustained  for  the  time  t^,  the  angular  displacement 
would  only  be  2t.ll.  This  has  no  significant  effect  on  the  analysis. 
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The  solution  for  9(t)  has  the  intuitively  expected 
characteristic  of  slow  early  growth  with  an  eventual  rapid  buildup  of 
acceleration.  For  the  example  in  Fig  B.l,  a/L  =  0.032,  the  following 
values  are  obtained: 


t/t„ 

Q(radians) 

0.127 

0.01 

0.293 

0.032 

0.508 

0.1 

0.824 

1.0 

2 

For  L  =  15.9mm  and  ou  =  10  kbars,  typical  of  conditions  nearer  the 
source,  t^  =  89  us  for  a  steel  helix  =  7.8)  and  t^  =32  us  for 
0u  =  1  (a  plastic).  This  may  explain  why  the  steel  helix  in  the 
first  300  mm  (LS  I V— 22 )  was  less  effective  than  the  polyolefin  helix 
in  the  first  229  mm  (LS  1 1 1—16) . 

2 

For  conditions  near  the  witness  plate  ou  =  3  kbars  is  more 
typical,  according  to  both  calculations  and  data.  The  plasma  leads 
the  particulates  by  about  60  uS.  t^  is  only  58  ys  for  =  1. 
therefore,  a  polyolefin  spiral  should  effectively  occlude  the  pipe. 

For  o,,  =  7.8,  t,  =  162  ys.  In  60  ys  a  steel  helix  with  a/L  *  0.032 
will  turn  only  2.8°,  allowing  the  edge  to  protrude  only  about  1.3  mm 
into  the  stream  (including  the  0.5  mm  spiral  thickness).  Although  the 
3.2  mm  lead  baffle  ( LS  II-5,  II 1-12 ,  and  IV-10,  11)  was  sufficient  to 
stop  the  particulates,  the  baffle  may  have  been  displaced  or  rotated 
into  the  stream;  hence,  a  1.3  mm  protrusion  of  spiral  may  be  marginal 
or  worse. 

Another  phenomenon,  not  yet  considered,  may  greatly  increase  the 
protrusion  of  the  helix.  A  twist  imparted  locally  will  propagate 
downstream  along  the  helix,  giving  it  an  initial  displacement,  a  head 
start.  If  the  head  start  is  only  0  =  a/L,  a  time  t  =  0.29  t^  can  be 
saved.  This  would  suggest  that  after  exposure  to  3  kbar  for  60  ys  in 
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the  above  example  a  steel  helix  might  be  turned  at  an  angle  appro¬ 
priate  for  t/t^  =  (60/162)  +  0.29  or  0.66,  i.e.,  14°,  and  protrude 
about  4  mm  into  the  stream.  The  propagation  of  twist  may  play  a 
significant  role  in  the  success  of  the  helix.  If  this  is  so,  the  use 
of  thin  rings  would  not  be  as  effective,  especially  if  the  material  is 
steel . 


It  remains  to  justify  the  neglect  of  the  strength  of  the  helix 
material.  The  sketch  in  Fig.  B.2  is  helpful  in  visualizing  the 
effect.  The  helix  is  viewed  as  uncoiled  into  a  straight  strip  but 
experiencing  a  twist  about  the  centerline  of  that  configuration.  The 
twist  considered  is  only  that  which  results  from  the  difference  between 
aerodynamical ly-induced  rotations  at  different  axial  locations.  The 
question  is  whether  the  torsion  resulting  from  this  twist  will  signifi¬ 
cantly  reduce  the  rotation.  The  twist  inherent  in  the  forming  of  the 
helix  is  ignored.  Bending  of  the  helix  in  the  plane  of  Fig.  8.1  is 
not  considered;  if  it  occurs,  it  only  enhances  the  intrusion  of  the 
helix  into  the  pipe. 

As  the  spiral  material  twists,  it  develops  internal  torsional 
stress  which  resists  further  twisting.  The  torsional  moment  about  the 
centerline  is 


T  = 


£ 

l 


where  E  =  shear  modulus  =  80  kbars  for  steel;  J  =  polar  moment  of 

s  3 

inertia  of  the  cross  section  about  the  centerline  =  aL  / 12 ; 

A  =  cross  sectional  area,  La.  Therefore, 


t  = 


B  3 

IT 


a3L  E. 


(B.8) 
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Consider  only  the  influence  of  the  torsion  resulting  from  the 
aerodynamic  twisting:  s  is,  then,  the  difference  in  8  over  a  length 
along  the  spiral.  A  generous  value  of  a /i  is  that  which  corresponds 


to  a  full  right  angle  twist  in  one 
e/£  =  cos  a/4R 

where  a  =  pitch  angle  and  R  =  pipe 

The  ratio  of  the  aerodynamic 
torsional  moment  per  unit  length  i' 


loop  around  the  pipe: 

(B.9) 

radius. 

torque  per  unit  length  to  the 
,  from  Eqs.  (B.2),  (B.8),  and  (B.9): 


i^>  <4  ( 

J  cs 


u2  /St"  9  cos  8 

3  ''  COS  a 

A 


Taking  0u2  =  10  kbar,  a  =  40°,  and  the  values  of  a,  L,  and  R  from 
the  PI  simulations. 


Ta  jlO5  for  6  =  ir/2 

^m  llO2  for  9  =  0 

The  higher  value  is  much  more  relevant  to  the  situation  posed.  In 
either  case,  however,  the  resistance  of  even  a  steel  spiral  to 
twisting  appears  to  be  small  compared  to  the  applied  aerodynamic 
moment,  and  the  approximations  made  in  the  dynamic  analysis  appear 
justified. 
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